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Abstract. The purpose of this paper is to describe a method for determining a cooling template 
for galaxies, using nearby galaxies, and applicable to future observations of distant galaxies. We 

observed two starburst galaxies (NGC 253 and Henize 2-10) with the Caltech Submillimeter Observatory in 
the rotational lines of carbon monoxide 12 CO(J=3-2), (J=6-5) and (J=7-6) for both, and also 12 CO(J=4-3) 
and 13 CO(J=3-2) for Henize 2-10 and in the 3 P2- 3 Pi fine-structure transitions of atomic carbon [CI] at 809 
GHz for NGC 253. Some of these observations have been made previously, but the present multitransition 
study (including data found in the literature) is the most complete to date for the two galaxies. From these 
observations, we have derived the properties of the warm and dense molecular gas in the galaxy nuclei. We 
used an LTE analysis and an LVG radiative transfer model to determine physical conditions of the interstellar 
medium in both sources and predicted integrated line properties of all CO transitions up to 12 CO(15-14). We 
£ — ■ found the observations to be in good agreement with a medium characterized by Tu ~ 50 — 100 K, v/qq ~ 30, 

n(H 2 ) > 10 4 cm" 3 and N( 12 CO) = 3.5 ± 1.0 x 10 18 cm" 2 for Henize 2-10 and characterized by T k « 70 - 150 K, 
^gg » 40, n(H 2 ) > 10 4 cm -3 and N{ 12 CO) = 1.5 ± 0.5 x 10 19 cm" 2 for NGC 253. A PDR model has also been 
used and here the data are well fitted (within 20 %) by a model cloud with a gas density of n(H)= 8. Oil. Ox 10° 
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cm" 3 and an incident FUV flux of X -20000 for Henize 2-10. For NGC 253, we deduced n(H)= 3.0±0.5x 10 5 
l' cm -3 and x ~20000 for the modelled cloud. The physical properties of warm gas and CO cooling curves 

1~ \ of the target galaxies are compared with those measured for the nucleus of the Milky Way and 

the Cloverleaf QSO. The gas properties and CO cooling curve are similar for the two starburst 
galaxies and the Cloverleaf QSO while the Milky Way nucleus exhibits lower excitation molecular 
03 '. g as - 
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1. Introduction to see how they relate to each other and to dust emis- 
sion (Gerin & Phillips 1998, 2000). However, to estimate 

Over the past decade, observations of distant galaxies have the relative contributions of the various species to the gas 

become possible in the millimeter and submillimeter bands coonri g ; it is often necessary to make assumptions for 

and considerable progress is expected for the next decade. ^he p Qwer m the various unobserved lines 
Much of the energy provided by star formation in those 

galaxies is down-converted to the submm and appears as It is well known that the fine-structure lines of ionized 
dust emission and the familiar gas cooling lines. Because carbon [CII] and atomic oxygen oxygen [01] contribute 
the information on distant galaxies is scarce, due to the most of the cooling of the neutral interstellar gas in galax- 
lack of resolution and sensitivity, it is essential to under- ies. These lines trace the cooling of the diffuse neutral gas 
stand nearby galaxies in order to provide templates for with a significant contribution from Photo-Dissociation 
more distant objects. In recent studies of emission of fine- Regions (PDRs). Fine-structure lines of ionized oxygen 
structure transitions of atomic carbon [CI] , ionized atomic [OIII] and nitrogen [Nil] trace the ionized gas, either in 
carbon [CII] and carbon monoxide (CO), we have begun HII regions([OIII]) or the diffuse ionised gas ([Nil]). In 
molecular gas, the cooling radiation is due to atomic car- 
Send offprint requests to: E.Bayet, e-mail: bon, carbon monoxide and water lines (although water 
estelle . bayet@lra . ens . f r (like most molecules) is not sufficiently widespread to con- 
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tribute much cooling to a galaxy as a whole). These pre- 
dictions have been confirmed by the COBE-FIRAS ob- 
servations of the Milky Way: apart from [CII] and [Nil] 
and probably [01], the most intense submillimeter lines 
are from C and CO (see Fixsen et al. 1999 and Sect. 14. 2[) . 
The relative contributions of the different lines of C and 
CO vary along the Galactic plane. Also C contributes less 
in proportion to the total cooling towards the Galactic 
Center than towards the rest of the disk (Bennett et al 
1994 and Fixsen et al. 1999). 

The detection of CO lines of distant galaxies remains 
an extremely difficult task and will continue so until 
ALMA is ready. However, detecting their submillimeter 
continuum emission is now feasible, with the advent of sen- 
sitive bolometer arrays. Due to the negative K-correction 
effect, the dust emission of distant galaxies piles up at sub- 
millimeter wavelengths irrespective of redshift. Additional 
information is needed to be able to estimate a photometric 
redshift, or to classify the galaxies. For instance, Yun & 
Carilli (2002) have proposed to use the cm radio contin- 
uum emission for obtaining photometric redshifts. From 
the model developed by Lagache, Dole & Puget (2003) 
to explain the evolution of number counts in the submil- 
limeter and far infrared (FIR), it appears that two main 
galaxy populations contribute to these number counts: a 
relatively local population of "cold" galaxies, similar to 
our own galaxy, in which the continuum emission is dom- 
inated by extended disk-type emission, and an additional, 
strongly evolving population of infrared (IR) luminous 
galaxies , similar to local starburst galaxies. From the in- 
formation gathered so far on high-J CO emission in galax- 
ies (this work; Israel & Baas 2002, 2003; Bradford et al. 
2003; Ward et al. 2003 and references therein), it is likely 
that high-J CO lines ( 12 CO(4-3) and above) will be more 
prominent in actively star forming galaxies compared to 
the less active galaxies in which only low-J CO lines and 
the carbon lines are expected. In this paper, from an 
extensive set of observations of submillimeter CO 
lines, we conclude that, relative to 12 CO(1-0), star- 
burst galaxies do exhibit stronger high-excitation 
CO lines ( 12 CO(4-3) and above), than the Milky 
Way galaxy. 

Determining the physical conditions of molecu- 
lar clouds in external galaxies can be difficult. For 
single-dish observations of external galaxies, the 
beam size is usually larger than the angular diam- 
eter of typical molecular clouds, hence the cloud 
emission is beam-diluted. In addition, molecular 
clouds are known to be clumpy (e.g. Stutzki &: 
Guesten 1990). The intensities of CO rotational 
lines are sensitive to the gas physical conditions 
(H2 density n(H2), and kinetic temperature T k ), 
the 12 CO column density N( 12 CO) as well as the 
velocity field. In this paper, we use several meth- 
ods for determining the 12 CO column density and 
the physical conditions, and for predicting the in- 
tensity of all CO rotational lines up to 12 CO(15- 
14). The LTE analysis is useful for constraining 



the kinetic temperature. Using LVG models, three 
parameters (n(H2),Tfc, N( 12 CO)) can be measured. 
PDR models include a self consistent treatment of 
the thermal and chemical processes and provide 
therefore a better theoretical understanding of the 
gas properties. 

We recall that, assuming LTE, the brightness temper- 
ature, Tg, of CO lines as a function of the kinetic temper- 
ature Tk can be written as: 

T B = FF {J u (T k ) - J v {T hg )) (1 - e- T ) (1) 

where FF is the source surface filling factor in the beam, r 
is the line opacity, Tb g is the CMB temperature {Tf, g = 2.7 
K) and J v (T k ) is described in Kutner & Ulich 1981 by: 

where fcf, is the Boltzmann's constant and h is the Planck 
constant. We show in Fig. ^ the brightness temperature, 
Tb, as a function of the kinetic temperature In this 
example, typical figures for Galactic molecular clouds have 
been chosen with FF — 1, N( 12 CO) = 1 x 10 18 cm" 2 and 
Aw = 1 ferns -1 , resulting in saturated CO lines. 

For a gas temperature around T k 10 K, low-J 
CO lines ( 12 CO(1-0), 12 CO(2-l) and 12 CO(3-2)) are eas- 
ily detected while mid and high-J CO lines ( 12 CO(5-4), 
12 CO (6-5) and above) need warmer gas to be detected 
(T fc > 15 K). 

Since we want to sample warm gas, we have observed, 
for both galaxies, the high-J CO lines: 12 CO(4-3) (only for 
Henize 2-10), 12 CO(6-5) and 12 CO(7-6) which are more 
sensitive to warm, dense gas directly involved in the star- 
burst. We also observed 12 CO(3-2) for both galaxies and 
13 CO (3-2) for Henize 2-10, to better constrain the models 
(see Sect. 15. l| l. 

Atomic carbon has also proved to be a good tracer of 
molecular gas (Israel & Baas 2001 and Gerin & Phillips 
2000) . The fine structure transition of atomic carbon 3 P2- 
3 Pi[CI] at 809 GHz for NGC 253 has also been observed 
here, since the intensity ratio between this line and the 
3 Pi- 3 Po[CI] at 492 GHz is a sensitive tracer of the total 
gas pressure. For NGC 253, Bradford et al. (2003) and 
Israel, White & Baas (1995) give the intensity of the fine 
structure transition of atomic carbon 3 Pi- 3 Po[CT] at 492 
GHz. For Henize 2-10, Gerin & Phillips (2000) give 
the intensity of the 3 P!- 3 P [CI] at 492 GHz but we 
could not find a measure of 3 P 2 - 3 Pi[CI] at 809 GHz 
to compute the line ratio. 

In the following, we present new CO and CI data for 
NGC 253 and for Henize 2-10 (see Sect.©. The galaxies 
are described in Sec. [21 the observations in Sec. The 
results and data analysis providing ISM properties 
are discussed in Sec. l4land in Sec. [SJ respectively, 
with conclusions in section 6. 
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Table 1. Basic properties of Henize 2-10 and NGC 253. 



T k (K) 

Fig. 1. Brightness temperature Tb (equations ^ 
and [5] in Sect. ^) vs. the kinetic temperature, 
Tk, assuming a unity filling factor (FF = 1). 
The first curve on the left is 12 CO(1-0) (full 
line) followed by 12 CO(2-l) (triangles), 12 CO(3-2) 
(squares), 12 CO(4-3) (full line)... up to 12 CO(15- 
14). For this plot, we used N(CO)= 1 x 10 18 cm -2 
and At; = 1 kms^ 1 resulting in r >> 1. 



2. The sample: 

Henize 2-10 is a blue compact dwarf galaxy with a metal- 
licity of 12+log(§) w 8.93 (Zaritsky, Kennicutt & Huchra 
1994). It presents a double core structure suggesting a 
galaxy merger (Baas, Israel & Koornneef 1994). From op- 
tical, infrared and radio observations of Henize 2-10, the 
derived distance is 6 Mpc (Johansson 1987). Henize 2- 
10 is rich in neutral atomic hydrogen (Allen, Wright & 
Goss 1976) and harbors a considerable population of Wolf- 
Rayet stars (D'Odorico, Rosa & Wampler 1983; Kawaraet 
al. 1987). Previous CO observations have been reported by 
Baas, Israel & Koornneef (1994); Kobulnicky et al. (1995 
and 1999) and Meier et al. (2001). The starburst is con- 
fined in the galaxy nucleus within a radius of 5" (150 pc). 
This galaxy was chosen because it shows intense and nar- 
row CO lines. 

NGC 253 is a highly inclined galaxy (i=78°, Pence 1980) of 
type Sc. With M82, it is the best nearby example of a nu- 
clear starburst (Rieke, Lebofsky & Walker 1988). Blecha 
(1986) observed 24 of the brightest globular clusters and 
estimated a distance of 2.4 to 3.4 Mpc, while Davidge 
& Pritchet (1990) analysed a color-magnitude diagram of 
stars in the halo of NGC 253 and conclude that the dis- 
tance is 1.7 to 2.6 Mpc. In the following, we shall use 
a value of D = 2.5 Mpc as in Mauersberger et al. (1996). 
NGC 253 has been extensively observed at all wavelengths 
(Turner 1985 (map in all four 18 cm lines of OH in the 
nucleus); Antonucci & Ulvestad 1988 (deep VLA maps of 
NGC 253 which show at least 35 compact radio sources); 
Carlstrom et al. 1990 (the J = 1-0 transitions of HCN 
and HCO + and the 3 mm continuum emission); Telesco 
& Harper 1980 (30 and 300 /jm emission); Strickland et 
al. 2004 (X-ray observations)). NGC 253 presents a very 
active and young starburst. It is among the first galaxies 
to have been detected in submillimeter lines and contin- 
uum. The NGC 253 nucleus has been mapped previously 
in various lines of CO and C (Bradford et al. 2003; Israel 



HENIZE 2-10 

RA(1950) 
DEC(1950) 
Distance 
Velocity (LSR) 
Optical size 
Position Angle 
Inclination 
Metallicity 

Absolute Luminosity Lb 



08 : 34 : 07.2 

-14 : 26 : 06.0 

QMpc 1 

+850 fcms" 1 

30" x 40" 

130° 

44° 

\2 + log{%) « 



.93 2 



1.7 x 10 ; 
1.2 x 10 ! 



L 
L 



NGC 253 



RA(1950) 
DEC(1950) 
Distance 
Velocity (LSR) 
Optical size 
Position Angle 
Inclination 
Metallicity 

Absolute Luminosity Lb 



00 : 45 : 05.7 

-25 : 33 : 38.0 

2.5Mpc 5 

+240 ferns -1 

27.5' x 6.8' 

51° 

78.5° 

12 + log(§) =8.99±0.31 2 



1.0 x 10 1U L 
2.7 x 10 10 L 



1 : Johansson 1987 
2 : Zaritsky et al. 1994 
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Computed from formula in Imanishi & Dudley (2000) and 
from Sanders & Mirabel (1996) 
4 : Baas, Israel & Koornneef 1994 

: As in Mauersberger et al. 1996 (see text) 
6 : Pence 1980 



& Baas 2002; Sorai et al. 2000; Harrison, Henkel & Russell 
1999; Israel, White & Baas 1995; Wall et al. 1991; Harris 
et al. 1991). The most intense emission is associated with 
the starburst nucleus with additional weak and extended 
emission from the disk (Sorai et al 2000). Both galaxies 
harbor superstar clusters in their nuclei (Keto et al. 1999) 
similar to the superstar clusters found in interacting galax- 
ies (Mirabel et al. 1998). Basic properties for both galaxies 
are summarized in Tabled 

3. Observations 

The observations were made during various sessions at 
the Caltech Submillimeter Observatory (CSO) in Hawaii 
(USA) with the Superconducting Tunnel Junction re- 
ceivers operated in double-side band mode. The atmo- 
spheric conditions varied from good (T225 < 0.1) to excel- 
lent (T225 ~ 0.06). We used a chopping secondary mirror 
with a throw of 1 to 3 arcmin, depending on the size of 
the source, and with a frequency around 1 Hz. Henize 2- 
10 is point-like for all lines. The NGC 253 nucleus has a 
size of 50" (see Fig. |3J) which is smaller than the chop- 
ping throw. For this galaxy, we used a 3' chopping throw 
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in 12 CO (3-2). There is no sign of contamination by emis- 
sion in the off beams. We restricted the chopping throw 
to 1' for 12 CO(6-5) and 12 CO(7-6) as the emission is very 
compact in these lines. Spectra were measured with two 
acousto-optic spectrometers (effective bandwidth of 1000 
MHz and 500 MHz). The first one has a spectral resolu- 
tion about 1.5 MHz and the second one about 2 MHz. 
The IF frequency of the CSO receivers is 1.5 GHz. The 
main beam efficiencies (77) of the CSO were 69.8%, 74.6%, 
51.5%, 28% and 28% at 230, 345, 460, 691 and 806 GHz 
respectively, as measured on planets. We used the ratio ~ 
to convert T^ into T mb . The beam size at 230, 345, 460, 
691 and 806 GHz is 30.5", 21.9", 14.5", 10.6" and 8.95" 
respectively 1 . 

The pointing was checked using planets (Jupiter, Mars 
and Saturn) and evolved stars (IRC 10216 and R-Hya). 
The pointing accuracy is around 5" . The overall cali- 
bration accuracy is 20%. Data have been reduced using 
the CLASS data analysis package. The spectra have been 
smoothed to a velocity resolution of 10 kms^ 1 and linear 
baselines have been removed. 

For NGC 253, we obtained spectra towards the central 
position for both 12 CO (7-6) line and the atomic carbon 
CI( 3 P2- 3 Pi) line simultaneously. We observed 70 positions 
in the 12 CO(3-2) line and 20 positions for the 12 CO(6-5) 
line. For Henize 2-10 we observed only the nucleus for the 
12 CO(3-2), 12 CO(4-3), 12 CO(6-5), 12 CO(7-6) and 13 CO(3- 
2) lines. The carbon line 3 Pi- 3 P [CI] is from Gerin & 
Phillips (2000). 



4. Data analysis 

4.1. Spectra and maps: 

Henize 2-10 and NGC 253 spectra at central positions are 
shown in Fig. Fig. shows 12 CO(3-2) and 12 CO(6-5) 
maps of NGC 253. We detected 12 CO(3-2), 12 CO(6-5) and 
12 CO(7-6) in both galaxies. For these positions, line inten- 
sities (A in Kkms -1 and I in Wm^sr^ 1 ) and line fluxes 
(F in WmT 2 ~) resulting from Gaussian fits are listed in 
Tableland in Table|3| together with the beam size of the 
observations. Because of likely pointing offsets (w 5" 
from pointing scans on planets and evolved stars) 
between different observing runs, we believe that 
differences in the peak positions of 12 CO(3-2) and 
12 CO(6-5) in NGC 253 are partly due to pointing 
errors. But we can not exclude that a difference 
in peak positions remains after the pointing cor- 
rection (see Fig. [3J). Three peaks appear on the 
high resolution 12 CO(1-0) maps obtained with the 
Nobeyama Interferometer (NRO), one near the 
center and two on each side at around 3 " south- 
west and 8 " north-east respectively (Paglione et 
al, 2004). The offset peak in the 12 CO (3-2) map 
may be associated with the north-east secondary 
peak. Tableland Table|3]also report relevant data found 
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Fig. 2. Spectra of Henize 2-10 (top six) and NGC 253 
(bottom four) at central positions as given in Tabled 
Velocities (horizontal axis) are given relative to 
the LSR (Vlsr) m kms^ 1 and the line intensities 
(vertical axis) are in units of T m & (K). The grey 
curves are Gaussian fits. For NGC 253, the line on 
the right hand side of CI( 3 P2- 3 Pi) is the edge of 
12 CO(7-6) line seen in the image side band. 



in the literature. For some spectra, a Gaussian may not 
represent the true line profile but given the weakness of 
the signal, fitting anything more complicated is not war- 
ranted. 

To be able to compare intensities of the different CO lines, 
we have also computed the line intensities (A and I) and 
fluxes (F) for a common beam size of 21.9". which is the 
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HENIZE 2-10 














Transition 


Freq 


Beam Size 


/ {Tmbdv) 


Intensity 


Flux 


References" 




(GHz) 


(") 


(Kkms" 1 ) 


(Wm- 2 sr- 1 ) 


(Wm" 2 ) 




CI( 3 Pi- a P ) 


492.162 


14.55 


4.2±0.8 


5.1xl0~ 10 


3.0xl0~ 18 


4& 5" 






21.90 


2.5±0.5 


3.0xl0~ 10 


3.8xl0 -18 


4 & 5 6 


i2 co(i-o) 


115.271 


40.00 


10.0±0.8 


1.6xl0 _il 


6.7xl0~ 19 


1 






21.90 


27.3±2.2 


4.3xl0 -11 


5.5xl0 -19 


1 






55.00 


4.9±0.2 


7.6xl0" 12 


6.1xl0" 19 


2 






21.90 


23.9±1.0 


3.7xl0 -11 


4.8xl0" 19 


2* 


12 CO(2-l) 


230.538 


21.00 


17.3±1.4 


2.2xl0~ 10 


2.5xl0 -18 


1* 






27.00 


6.8±0.8 


8.5xl0 -11 


1.7xl0" 18 


2 






21.90 


9.4±1.1 


1.2xl0~ 10 


1.5xl0" 18 


2 


12 CO(3-2) 


345.796 


21.90 


11.5±2.3 


4.9xl0~ 10 


6.2xl0~ 18 


5* 






21.00 


23.2±2.1 


9.8xl0" 10 


1.2xl0" 17 


1 






22.00 


16.6±0.6 


7.0xl0" 10 


9.0xl0" 18 


3 


12 CO(4-3) 


461.041 


14.55 


18.6±4.2 


1.9xl0~ 9 


l.lxlO" 17 


5 






21.90 


10.9±2.4 


l.lxlO" 9 


1.4xl0~ 17 


5* 


12 CO(6-5) 


691.473 


10.60 


15.7±3.1 


5.3xl0~ 9 


1.6xl0~ 17 


5 






21.90 


6.8±1.3 


2.3xl0~ 9 


2.9xl0~ 17 


5* 


12 CO(7-6) 


806.652 


8.95 


15.2±3.0 


8.2xl0" 9 


1.7xl0~ 17 


5 






21.90 


5.8±1.2 


3.1xl0 -9 


4.0xl0~ 17 


5* 


"CO (1-0) 


110.201 


40.00 


< 0.5 


< 6.8x10"" 


< 2.9xl0 _aj 


1 






21.90 


< 1.4 


< 1.8xl0~ 12 


< 2.4xl0~ 20 


1 






57.00 


0.3±0.1 


3.9xl0~ 13 


3.4xl0~ 20 


2 






21.90 


1.6±0.5 


2.2xl0~ 12 


2.8xl0" 20 


2* 


13 CO(2-l) 


220.399 


21.00 


0.9±0.2 


9.9xl0~ 12 


1.2xl0" 19 


1* 


13 CO(3-2) 


330.588 


14.00 


2.3±0.6 


8.5X10" 11 


4.4xl0 -19 


1 






21.90 


1.3±0.3 


4.8xl0~ n 


6.1xl0 -19 


1* 






21.90 


1.8±0.5 


6.7xl0~ n 


8.5xl0" 19 


5 



a References: 1: Baas, Israel & Koornneef 1994; 2: Kobulnicky et al. 1995; 3: Meier et al. 2001; 4: Gerin & 
Phillips 2000; 5: this work; *: used for models. b We analysed again spectra from Gerin & Phillips 2000. 



size of the 12 CO(3-2) beam at the CSO. They are listed 
in Tableland in Table |3 and identified with an asterisk 
in column 7 (named References). To do so, we have used 
the following assumptions: 

— Henize 2-10: we modelled the source with an axisym- 
mctric Gaussian distribution, with a full width half 
maximum (FWHM) of 13", as found by Meier et al. 
2001. 

— NGC 253: we made use of the 12 CO(6-5) map to de- 
termine the source size at 690 GHz and model the 
emission at higher frequencies. We modelled the nu- 
cleus as an elliptical source with Gaussian intensity 
profiles. The 12 CO (6-5) map can be well fitted with 
such a model, with half maximum width of 23" along 
the major axis and 11" along the minor axis, which is 
very similar to the size of the CS(J=2-1) emission in 
the nucleus (Peng et al. 1996). The source appears to 
be resolved along the major axis but unresolved along 
the minor axis. We used the same source model for 
12 CO(7-6) and CI( 3 P 2 - 3 Pi), as it fits our restricted 
data set. 

— I is derived using the formula: 

2 x kb x v 3 



I = 



x A 



(3) 



liWm^sr- 1 ) = 1.02xl0" 18 x(-^— ) 



3 ( 

X 



A 



\Kkms~ 



(4) 



where c is the speed of light, v is the line frequency in 
GHz and A is the line area in Kkms^ 1 . 
To derive the flux, F in Wm" 2 , we multiply I in 
Wm~ 2 sr~ 1 by the beam solid angle fl(B) in sr de- 
fined as: 



n(B) = 1.133 x B 2 x 



1 



206265 2 



(5) 



where B is the half power beam width (HPBW) in 
arcsec. The estimated error for our data in Table [21 
and in Table 01 is about 20%. 



4.2. C and CO cooling: 

Our observations are designed to provide essential 
information on the cooling and consequently on the 
thermal balance of the interstellar medium in Henize 2-10 
and NGC 253. We shall also deduce which CO line(s) 
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Table 3. Results for the nucleus of NGC 253. 
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Freq 


Beam Size 


/ (Tmbdv) 
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(GHz) 
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(Kkms- 1 ) 


(Wm- 2 sr- 1 ) 


(Wm" 2 ) 




CI( 3 Pi- 3 P ) 


492.162 


10.20 


575.0±115.0 b 


7.0xl0~ 8 


1.9xl0~ 16 


5 






21.90 


204.0±41.0 


2.5xl0" 8 


3.2xl0" 16 
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43.00 


98.0±19.6 b 
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5.9xl0 -16 


5 






22.00 


290.0±45 


3.5xl0" 8 


4.6xl0" 16 


10* 






23.00 


320.0±64.0 6 


3.9xl0" 8 


5.5xl0" 16 


4 


CI( 3 P 2 - 3 Pi) 


809.902 


8.95 


188.5±37.7 


1.0xl0~ 7 


2.2xl0~ 16 


11 






21.9 


58.1±11.6 


3.2xl0~ 8 


4.0xl0~ 16 


11* 


12 CO(1-0) 


115.271 


43.00 


343.0±68.6 b 


5.4xlO~ 10 


2.6xl0~ 17 


2 






23.00 


920.0±82.8 


1.4xl0" 9 


2.0xl0~ 17 


6* 


12 CO(2-l) 


230.538 


23.00 


1062.0±116.8 


1.3xl0~ 8 


1.9xl0~ 16 


7* 






21.00 


926.0±185.2 i ' 


1.2xl0" 8 


1.4xl0~ 16 


2 


12 CO(3-2) 


345.796 


21.90 


815.6±163.1 


3.5xl0 -8 


4.4xl0 -16 


11* 






23.00 


998.0±139.7 


4.2xl0" 8 


5.9xl0~ 16 


7 






23.00 


1194.0±238.8 6 


5.0xl0 -8 


7.1xl0 -16 


2 






14.00 


1200.0±240.0 b 


5.1xl0" 8 


2.6xl0" 16 


5 






21.90 


618.9±123.8 


2.6xl0~ 8 


3.3xl0~ 16 


5 






22.00 


680.0±60.0 


2.9xl0" 8 


3.7xl0" 16 


8 


12 CO(4-3) 


461.041 


15.00 


507.0±101.4 i ' 


5.1xl0~ 8 


3.0xl0~ 16 


3 






21.90 


285.3±57.1 


2.9xl0~ 8 


3.7xl0 -16 


3 






22.00 


1019.0±120 


l.OxlO" 7 


1.3xl0" 15 


9* 






10.40 


2160.0±432.0 6 


2.1xl0~ 7 


6.2xl0 -16 


5 






21.90 


786.2±157.2 


7.9xl0~ 7 


l.OxlO -15 


5 


12 CO(6-5) 


691.473 


10.60 


1394±278.8 


4.7xl0" 7 


1.4xl0" 15 


11 






21.90 


518.3±104.3 


1.8xl0" 7 


2.2xl0" 15 


11* 






8/30 


861±258.3 


2.9xl0" 7 


5.0xl0" 15 


1 


12 CO(7-6) 


806.652 


8.95 


810.2±162.0 


4.3xl0 -7 


9.3xl0 -16 


11 






21.90 


249.9±50.3 


1.3xl0" 7 


1.7xl0~ 15 


11* 






11.5/60 


1370±411 


7.3xl0 -7 


2.6xl0~ 15 


10 


i3 CO(l-0) 


110.201 


23.00 


80±8.0 


l.lxl0~ iu 


1.5xl0~ 18 


7* 


13 CO(2-l) 


220.399 


23.00 


82±9.8 


9.0xl0" 10 


1.3xl0~ 17 


7* 






21.00 


104±20.8 i ' 


l.lxlO" 9 


1.3xl0~ 17 


2 


13 CO(3-2) 


330.588 


23.00 


90±12.6 


3.3xl0 -9 


4.7xl0~ 17 


7* 






23.00 


210±42.0 b 


7.8xl0" 9 


l.lxlO -16 


2 



a References: 1: Harris et al. 1991; 2: Wall et al. 1991; 3: Guesten et al. 1993; 4: Harrison et al. 1995; 5: Israel, 
White & Baas 1995; 6: Mauersberger et al. 1996; 7: Harrison, Henkel & Russell 1999; 8: Dumke et al. 2001; 9: 
Israel & Baas 2002; 10: Bradford et al. 2003; 11: this work; *: used for models. 
6 The error is not given in the referenced paper. It has been estimated to be 20%. 



contribute the most to the total observed CO cooling 
and estimate the total observed cooling of C and CO by 
summing intensities in Wm~ 2 sr~ 1 of all transitions listed 
in Table and in Table EH with asterisks (both literature 
data and our dataset). We have computed the observed 
C and CO cooling in the galaxy nuclei for a beam size of 
21.9", this corresponds to linear scales of 640 pc and 270 
pc for Henize 2-10 and NGC 253 respectively. 
For Henize 2-10 and NGC 253 we measured a total 
observed CO cooling rate of 7.2 x 10~ 9 Wmr 2 sr~ x and 
4.6 x 10~ 7 Wm^sr' 1 respectively. For NGC 253 the 
lines contributing the most to the observed CO cooling 
are 12 CO(6-5) (39.2% of the total intensity) followed 
by 12 CO(7-6) (28.3%). For Henize 2-10, it is 12 CO(7-6) 
(43.1%) followed by 12 CO(6-5) (31.9%). 
The observed total cooling rates of neutral carbon 



C for Henize 2-10 and NGC 253 are respectively 
3.0 x lO" 10 WmT 2 sr- x (CI( 3 Pi- 3 P ) transition only) 
and 6.7 x 10~ 8 Wm^sr- 1 (for NGC 253, CI( 3 P 2 - 3 Pi) 
represents 52.2% of the total). 

These results show that 12 CO(6-5) and 12 CO(7-6) are 
contributing the most to the total observed CO cooling, 
with very similar percentages for both galaxies. It is 
natural to wonder whether 12 CO(5-4) and 12 CO(8-7) 
are strong also. This will be investigated in Sec. 15.21 and 
Sec.O 

In starburst nuclei, CO cooling is larger than C cooling by 
a factor of > 10, which explains why it is easier to detect 
CO than C in distant galaxies. Similar results have been 
obtained on J1148+5251 (z=6.42) and on PSS2322+1944 
(^=4.12) (Bertoldi et al (2003), Walter et al. (2003) 
and Cox et al. (2002)). Barvainis et al (1997) detected 
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CI( 3 Pi- 3 P ) (3.6 ± 0.4 Jykms" 1 ) in the Cloverleaf 
quasar at redshift z=2.5. In the same quasar, Wei/3 et al 
(2003) detected CI( 3 P 2 - 3 Pi) (5.2 ± 0.3 Jykms' 1 ) and 
12 CO(3-2) (13.2 ± 0.2 Jykms- 1 ). In distant objects, C 
cooling also seems to be weaker than CO cooling. 



5. CO models: 

In this section, we use the measured CO line ratios and 
intensities (I and A) to determine the physical conditions 
of molecular gas, namely the kinetic temperature, the 
gas density, the CO column density and the FUV flux: 
X- In the first section, we start with an LTE analysis. In 
Sect. 15.21 we use an LVG radiative transfer model and 
in Sect. 15.31 we discuss the use of a PDR model. In the 
last section, we discuss the similarities and differences 
between the two galaxies, the nucleus of the Milky 
Way and the distant QSO "the Cloverleaf". The 
latter has been chosen as being representative of 
distant, actively star forming galaxies. 



o - 



V 



nbeam 
C0(6-5) 



10 -10 
Minor Axis offset (arcsec) 



Fig. 3. Map of the velocity integrated intensity of 12 CO(3- 
2) (Top) and 12 CO(6-5) (Bottom) towards NGC 253 nu- 
cleus. The CO emission is integrated over the full velocity 
range - 100 to 700 fcms -1 and 60 to 400 fcms -1 , respec- 
tively. For both maps, missing grid positions are interpo- 
lated and the crosses show the observed positions. The 
first and the second contours of the 12 CO(3-2) map are of 
41.6 Kkms -1 and 81.6 Kkms" 1 respectively; others con- 
tours range from J T mb dv = 81.6 to 816 Kkms- 1 in steps 
of 81.6 Kkms -1 . For the 12 CO(6-5) map, contours are of 
43 Kkms' 1 , 1/11 z^'— - 1 



141 Kkr 



282 Kkms- 



564 Kkms , 

846 Kkms- 1 , 1128 Kkms- 1 and 1410 Kkms- 1 . The peak 
value of the 12 CO(3-2) map is / T mb dv « 816 Kkms- 1 . 
The peak value of the 12 CO (6-5) map is J T mb dv « 1400 
Kkms -1 . The maps are rotated by 51° to match the ma- 
jor axis position angle of the CO emission in NGC 253. 
The 12 CO (3-2) map has been shifted by -5" in mi- 
nor axis offset (see text in Sect.|3|and in Sect. HTT|) . 
The black box in the upper figure represents the 
size of the lower figure. 



5.1. LTE analysis: 

In the local thermodynamical equilibrium (LTE) approxi- 
mation, we assume that CO is thermalized, hence the rel- 
ative populations of its energy levels are functions of the 
kinetic temperature (assumed uniform) only. We discuss 
two limiting cases, applied to 12 CO and 13 CO: 

o the lines are optically thick (optical depth: t>1); the 
line intensity ratio (eg i2 CO (2~i) > can ^ e wr itten: 



R 



thick 
32/21 



FF x T 12 CO(3 - 2) _ J„ 32 (T fc ) - J^ 2 {T bg ) 
FF x T 12 CO(2 - 1) 



Jv2l(Tk) — Jv2l(Tbg) 

(6) 



-"-32/21 ~ 1 Ana -"76/65 
TDthick -i OT1 j jjthick 

-"65/32 ~ 1 ancl -"76/32 



with J v {T k ) defined as in Eq in Sec. □ In LTE, 

1 for T k > 30 K, while 

1 for T k > 100 K. Plots 
of the line intensity ratios (R thlck ^ as a function of 
the kinetic temperature Tk for the optically thick case 
in the LTE approximation are given in Fig. From 
this figure, we conclude that the line ratios combining 
high-J ( 12 CO(6-5) or 12 CO(7-6)) with low-J ( 12 CO(3- 
2)) lines are the most useful for constraining the kinetic 
temperature. 

o One line is optically thick and the other is optically 
thin; this is likely when comparing line intensities for 
different isotopologues of CO in the same J transition 
(eg 12 CO(3-2)/ 13 CO(3-2)). Assuming equal excita- 
tion temperatures for both species, optically thin 13 CO 
lines and optically thick 12 CO lines, we get: 



-^12/13 — 



1 



FF x T 12 CO _ A( 12 CO) ^ 
FF x T 13 CO ~ A{ 13 CO) ~ t( 13 CO) 



(7) 
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Fig. 4. Integrated intensity lines ratios, ii*"* , (denned as 
in Eq. EJ), vs kinetic temperature in K. We assume LTE, 
and optically thick lines (optical depth, t>1). 



X galaxy 

t( 12 CO) 



(8) 



where X ga i axy is the 12 CO/ 13 CO abundance ratio. 
Previous work, using LVG models, found, for Henize 
2-10, X ~ 30 (Baas, Israel k Koornneef 1994) and 
for NGC 253, X ~ 40 (Henkel & Mauersberger 1993; 
Henkel et al. 1993 and Israel & Baas 2002). From the 
literature data, summarized in TableEl and in Tabled 
we deduced t( 12 CO(1 - 0)) ~ 2.0 and r( 12 CO(2 - 
1)) ~ 1.6 for Henize 2-10. For NGC 253, we obtained 
r( 12 CO(l - 0)) ~ 3.5 and r( 12 CO(2 - 1)) ~ 3.1 (r 
increase for larger X ga i axy : assuming X NGC 253 — 50 
as in Bradford et al. (2003) r( 12 CO(l - 0)) ~ 4.3 and 
r( 12 CO(2-l))~3.9. 

Unfortunately, the opacity of the 12 CO(3-2) line for 
both galaxies is not accurately determined due to sys- 
tematic differences between data sets (Baas, Israel 
& Koornneef 1994; Harrison, Henkel & Russell 1999 
and this work). For Henize 2-10, r( 12 CO(3 - 2)) 
ranges between 1.7 and 3.4, and for NGC 253, it 
ranges between 3.0 and 5.8 for the abundance ra- 
tio given above (XnqC253 = 40). We note that the 
largest error is for the 13 CO integrated intensity which 
is the most difficult line to observe: for NGC 253, 
A( 13 CO(3-2)) varies between 90 Kkms' 1 (Harrison, 
Henkel & Russell 1999) and 210 Kkms' 1 (Wall et al. 
1991) while A( 12 CO(3-2)) varies between 619 Kkms' 1 
(Israel, White & Baas 1995) and 1194 Kkms- 1 (Wall 
et al. 1991) (at the same resolution). 

Because the observed ratio is -Rg5/32 ~ 0-6 for both Henize 
2-10 and NGC 253 (see Table 0}, we conclude, from the 
LTE analysis, that the 12 CO(3-2) line opacities are 
moderate (1 < r < 10) and the gas is warm (Tk > 30 K). 
To go further and constrain the gas density, we now use 
LVG models. 



5.2. LVG model: 

The detection of the high-excitation CO lines directly in- 
dicates the presence of large amounts of warm and dense 
gas in these galaxies. The CO excitation is modelled in 
terms of standard, one-component, spherical LVG radia- 
tive transfer model with uniform kinetic temperature and 
density (Goldreich & Kwan 1974; De Jong, Dalgarno & 
Chu 1975). Though the assumption of uniform phys- 
ical conditions is crude, these models represent a 
useful step in predicting intensities of the submil- 
limeter CO lines, since, in any case, we lack in- 
formation concerning the gas distribution in the 
galaxy, particularly its small scale structure. 

There are four main variables in LVG models: the CO 
column density divided by the line width: N( 12 CO)/Ai;, 
the molecular hydrogen density: n(H2), the kinetic tem- 
perature: Tfc and the vaqq abundance ratio: X ga iaxy The 
linewidth, Av, is constrained by the fit of the 
12 CO(3-2) spectra (FWHM). We have used Aw = 
60 ferns" 1 for Henize 2-10 and Av = 190 ferns" 1 for NGC 
253 (see Fig. 01. The ts§§ abundance ratio varies between 
30 and 40 for Henize 2-10 (X Henize2 -io) and between 30 
and 50 for NGC 253 {X NGC253 ). 

To make quantitative estimates of the physical 
parameters of the molecular gas in nuclei, we 
have compiled information about line emission for 
all CO rotational transitions observed so far and 
scaled them to a common beam size of 21.9" (see 
Table El and Table El and Sec. l4TT|) . 

Model solutions for the physical parameters, for both 
sources, are presented in Table Because acceptable fits 
can be obtained over a large domain of parameter space, 
we have chosen to present two models which bracket the 
range of possible solutions for the kinetic temperature. 
The first one has a "low kinetic temperature" (Tk < 
80 K), shown with crosses in Figs. El and EJ) and the 
second one has a "high Tk" (Tk > 100 K), shown with 
triangles in Figs. and E|) (see Table 0J. The gas den- 
sity is not well constrained but must be at least 10 cm~ 3 . 
The best fit 12 CO/ 13 CO value is 30 for Henize 2-10 and 40 
for NGC 253, both compatible with previous analyses. We 

2 to 1.0 x 10 20 cm- 2 



16 



cm 



cm 2 to 3.5 x 10 19 cm 



varied N( 12 CO) from 1.0 x 10 
for Henize 2-10 and 3.5 x 10 17 
for NGC 253. CO column densities are mainly constrained 
by the 13 CO data (see Table Hand Table 0). Acceptable 
fits are found : N( 12 CO)=3.5 ± 1 x 10 18 cm" 2 for Henize 



2-10 and for N( 12 CO) = 1.5 ± 0.5 x 10 



if) 



cm 



" 2 for NGC 



253. We determined the "best" fit using a least square 
fit method for the line area ratios listed in Table E] As 
we wanted to derive the properties of warm gas, 
we have given more weight to line ratios which in- 
clude high- J CO lines, e.g. isco^-l] ' ^cop-e) and 
Z C C of 7 ~-% (see Sect.QLU). 

Predicted line integrated areas, A in Kkms -1 and line 
intensities I in Wm^sr' 1 for all CO transitions to J=15- 
14, are also shown in Figs. Eland EJ The predicted line area 
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A is obtained by multiplying the antenna temperature, 
given by models, by the line width Av and by the surface 
filling factor, FF. For each observed CO line, we estimated 
the surface filling factor (FF) of molecular clouds in the 
beam using the ratio T mb (observations)/T mb (model). We 
made an average, weighted by the S/N, for the FF 
pertaining to each observed transition. For Henize 
2-10 and NGC 253, we obtained filling factors of FF = 
8.8 ±3 x 1(T 3 and fF = 8.7i3x 10" 2 respectively. The 
line intensities, I, are derived from the model A values 
through Eq.[3]and Eq.^J We give in Tablc^J line intensity 
ratios of CO transitions derived from observations and 
from models. 

The low-J CO transitions ( 12 CO(1-0) and 
12 CO(2-l)) are not well fit by the adopted LVG 
models, because we have chosen to constrain the 
models with high-J transitions of CO as we are in- 
terested in the warm gas properties. In order to 
study the properties of the low excitation molec- 
ular gas, we could have introduced another gas 
component in the models as in Harrison, Henkel 
& Russel (1999) or Bradford et al. (2003) for NGC 
253. Here, the LVG models predict too low an in- 
tensity for the 12 CO(1-0) and 12 CO(2-l) lines in 
this source, supporting the need for a low excita- 
tion gas component. For NGC 253, we obtained a 
similar kinetic temperature Tjc as in Bradford et 
al. 2003, but a lower density n{H2) and a higher CO 
column density N( 12 CO). For NGC 253, Paglione 
et al. (2004) found good fits with LVG models with 
N( 12 CO)/Aw between 4.0 xlO 16 cm^km^s and 1.3 x 
10 18 cm km s, in good agreement with our value 
(N( 12 CO)/Au = 8.0 x 10 16 cm- 2 km- 1 s). Recent mea- 
surements of rotational lines of H2 (Rigopoulou et al. 
2002) suggest that warm gas with Tk — 150 K is present 
in NGC 253 which is compatible with our models. 

The LVG models have been used for predict- 
ing the CO line intensities for 12 CO(1-0) up to 
12 CO(15-14). We have derived the total CO cooling 
from these predictions. For Henize 2-10, we obtain with 
the "low IV' model (T K = 50 K), 1.3 xlO" 8 Wm- 2 sr~\ 
and with the "high T K " model (T K = 100 K), 2.0 x 
10~ 8 Wm^sr- 1 . With the "low TV model, the lines 
which contribute the most to the total CO cooling are 
12 CO(6-5) (23.4%), followed by 12 CO(5-4) (19.7%) and 
12 CO(7-6) (19.5%). With the "high TV' model, the lines 
which contribute the most to the total CO cooling are 
12 CO(7-6) (19.1%) followed by 12 CO(6-5) (18.2%) and 
12 CO(8-7) (15.4%). 

For NGC 253, with the "low TV' model {T K = 70 K) 
we deduced a total CO cooling of 1.0 x 10~ 6 W mr" 1 sr~ x 
and the most important line is 12 CO(7-6) (21.1%) fol- 
lowed by 12 CO(6-5) (20.4%), 12 CO(8-7) (15.7%) and 
12 CO(5-4) (15.2%). With the "high TV model (T K = 
150 K), we obtained a total predicted CO cooling of 
1.5 x 10 -6 Wm~ 2 sr^ 1 , with the most intense lines be- 
ing 12 CO(8-7) and 12 CO(7-6) (16.4% for each) followed 
by 12 CO(6-5) (13.6%) and 12 CO(9-8) (13.4%). 



Table 4. Observed and predicted ("low TV and "high 
Tk" LVG models) line area ratios for Henize 2-10 and 
NGC 253. 



A 


NGC 253 


HENIZE 2-10 


[i\ KfflS ) 


observat ions 


observations 


LA one 3 .— 91 

12 CO{6-5) 


1.6±0.6* 


1.7±0.7* 


u CO(3-2) 
u CO(7-6) 


3.3±1.3* 


2.0±0.8* 


12 CO(6-5) 
u CO(7-6) 


2.1±0.8* 


1.2±0.5* 


12 CO(1-0) 

13 co(i-o) 


11.5±2.2* 


14.9±5.8* 


12 CO(2-l) 
13 CO(2-l) 


13.0±3.0* 


19.2±5.8* 


u CO(3-2) 
13 CO(3-2) 


9.1±3.1* 


8.9±3.8* 




"low T x " model 


"low T K " model 




1 K = < U il 


±K — t)U K 


LZ CO(3— 2) 
12 CO(6-5) 


1.6 


2.0 


u CO(3-2) 
12 CO(7-6) 


2.5 


3.8 


12 CO(6-5) 
12 CO(7-6) 


1.5 


1.9 


12 CO(1-0) 


25.5 


19.1 


12 CO(2-l) 
1*00(2-1) 


13.4 


9.9 


12 CO(3-2) 
13 CO(3-2) 


10.1 


8.1 




"high Tk" model 


"high T K " model 




Tk = 150 K 


Ik = 100 K 


12 CO(3-2) 
12 CO(6-5) 


1.4 


1.6 


12 CO(3-2) 
u CO(7-6) 


1.8 


2.5 


12 CO(6-5) 
u CO(7-6) 


1.3 


1.5 


12 CO(1-0) 
la CO(l-0) 


15.1 


23.9 


12 CO(2-l) 
13 CO(2-l) 


20.7 


14.8 


12 CO(3-2) 
13 CO(3-2) 


14.8 


9.9 



* : ratio derived from observations marked with asterisks in 
Tableland in Table 



12 CO(6-5) and 12 CO(7-6) therefore appear to con- 
tribute significantly to the CO cooling. Also, CO lines with 
Jupper ^ 10 are predicted to be weak and will not have 
significant antenna temperatures (see Figs. Eland El plots 
on the left side). In addition to 12 CO(6-5) and 12 CO(7-6), 
data for 12 CO (8-7) and 12 CO (9-8) would be most useful in 
discriminating between models, and for a more accurate 
determination of the CO cooling. 13 CO (6-5) data would 
also be extremely useful for constraining the models and 
for measuring the opacity of 12 CO(6-5) line. 

5.3. PDR model: 

To progress further in the analysis of the physical condi- 
tions in the starburst nuclei, we made use of PDR mod- 
els. Such models have been used here to under- 
stand the properties of the interstellar medium, 
since they take account of all the relevant physical 
and chemical processes for thermal balance. The 
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Fig. 5. LVG and PDR model calculations compared with observations for Henize 2-10. On the left side, we plot 



integrated intensities in Kkms vs J 



upper 



On the right side, we plot I in Wm 



vs J 



upper 



In both figures, 



grey crosses represent the "low temperature" LVG model while grey triangles represent the "high temperature" LVG 
model. The starred grey triangles represent results of the PDR model (see Sect. 15.3]) . Observations (with error bars) 
taken from literature and from our data set are shown in black. We have reported two observations for 12 CO(3 — 2) 
in Henize 2-10: the black cross is from our data and the black circle is from Meier et al. 1991. The last one seems to 
be agree better with the model predictions. As Henize 2-10 is a point source, a small pointing error can explain this 
difference. 



Table 5. Parameters of the "best" LVG models. 





Henize 2-10 


NGC 253 


n(H 2 ) (cm" 3 ) 


>1 x 10 4 


>1 x 10 4 


N( i2 CO) (cm"") 


3.5±1 x 10 1S 


1.5±0.5 x 10 ly 


vj c 


30 


40 


Av* (kmiT 1 ) 


60 


187 


Tk (K) 


50-100 


70-150 


FF 


8.8±3 x lO-* 


8.7±3 x 10-* 



* deduced from Gaussian fits to the 12 CO(3 — 2) line profiles. 

PDR models accurately reproduce the steep ki- 
netic temperature gradient near cloud edges, when 
illuminated by intense FUV radiation. Such models 
have been developed during the past two decades, for a va- 
riety of astrophysical sources, from giant molecular clouds 
illuminated by the interstellar radiation field to the condi- 
tions experienced by circumstellar disks, very close to hot 
massive stars (Tielens & Hollenbach 1985; Van Dishoeck 



& Black 1986, 1988; Wolfire, Hollenbach & Tielens 1990; 
Hollenbach, Takahashi & Tielens 1991; Abgrall et al. 1992; 
Le Bourlot at al. 1993; Koster et al. 1994; Sternberg & 
Dalgarno 1995; Draine & Bertoldi 1996, Stoerzer et al. 
1996, Lhuman et al. 1997; Pak et al. 1998; Hollenbach & 
Tielens 1999 and Kaufman et al. 1999). 

We adopted here the PDR model developed by 
Le Bourlot et al. (1993) for Galactic sources (see 
also Le Petit, Roueff & Le Bourlot 2002) . The 

source is modelled as a plane-parallel slab, illuminated 
on both sides by FUV radiation to better reproduce the 
starburst environment as massive stars and giant molec- 
ular clouds are spatially correlated. Model parameters 
include the gas density, assumed uniform, the intensity of 
the illuminating FUV radiation, the gas phase elemental 
abundances, the grain properties, the gas to dust ratio, 
etc. Because the metallicities of both NGC 253 and 
Henize 2-10 (see Table are close to solar, we used a 
model with Milky Way abundances (12 + log(%) =8.90 
± 0.04 Boselli, Lequeux & Gavazzi 2002). We have 
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LVG model: crosses 
Tk = 70 K 

N(C0) = 1.5*10 cm 2 - 

1500 - n(H2) = 1*10 4 cm" 3 _ 
LVG model: triangles 
" Tk = 150 K 

N(C0) = 1.5*10 19 cm 2 

n(H2) = 1*10 4 cm -3 



„1000 -t 



> 



a 500 
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HI 
Pi 



10" 



10" 



10 



10" 



10" 



10 



: PDR model: others symbols 
n(H) = 3*10 5 cm" 3 
X=20000 



NGC 253 



observations 
models 



5 10 15 

J(upper) 



Fig. 6. LVG and PDR models calculations compared with observations for NGC 253. The observations are shown in 
black with error bars while model predictions are shown in grey. See caption of Fig. 



adopted standard grain properties and gas to dust 
ratio appropriate for Galactic interstellar clouds. 

We have sampled a wide range of the parameter space, 
varying the gas density, n(H) and the incident FUV 
flux, x x Go, where Go is the local average interstellar 
radiation field (ISRF) determined by Draine (1978) 
(Go = 2.7 x 1CT 3 ergem^s- 1 ). The 12 C/ 13 C ratios 
(X He mze2-w and X NGC 253) are the same as for the 
"best" LVG models (see Table |3J). Models are constrained 
by the ratio of integrated intensities in ergcmT 2 s -1 sr -1 
listed in Tableland in Table {E{ergcmr 2 s~ x sr~ v ) = 
l{Wm- 2 sr- x ) x 10 3 ). 

Parameter pairs (n(H); x) which fit our dataset, 
are listed in Table El X seems to be constrained 
by the 12 CO(J+l^J)/ 13 CO(J+l-^J)(e.g. 12 CO(3- 
2)/ 13 CO(3-2)) ratios while n(H) is more sensi- 
tive to the line ratios involving two 12 CO lines: 
12 CO(J+l^J)/ 12 CO(J'+l^J')(e.g. 12 CO(3-2)/ 12 CO(6- 
5)) ratio. The difference between observations and model 
outputs is estimated to be of the order of 20%. Emissivity 
ratios obtained from the observations and from the 
models are listed in Table [7| 



Model predictions for the CO line emissivities are 
shown in Fig. and Fig. for Henize 2-10 and NGC 253 
respectively. The model predictions have been scaled to 
match the observed line intensities. As stated above, 
PDR models have been developed for local in- 
terstellar clouds. The velocity dispersion of the 
modelled cloud is a parameter in those models, 
which is used for computing the photo-dissociation 
rates of H2 and CO, and the line emissivities. 
This parameter is set to 1 ferns -1 , a typical figure 
for local molecular clouds (see Wolfire, Tielens & 
Hollenbach 1990). Using a PDR model for fitting 
the galaxy observations is complicated by the fact 
that, in a galaxy, many PDRs contribute to the 
signal detected in each beam, resulting in a broad 
line (tens to hundred of kms *), compared to a sin- 
gle PDR line (1 kms *). To correct for this effect, 
the model line emissivities have been multiplied by 
the line width ratio, , where Av(galaxy) as 

190 kms- 1 and 60 ferns" 1 for NGC 253 and Henize 
2-10 respectively, and Av(PDR) = 1 ferns -1 . Once 
this correction is performed, PDR model results 
are compared with observed data in the same way 
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as the LVG models (see paragraph 5 in Sect. IBT2*|) . 
and the surface filling factor of the emission in 
the beam, PDR_FF, is computed. The PDR_FF 
is 9.4x10 2 for NGC 253 and 1.5 xlO 3 for Henize 
2-10. For NGC 253, the surface filling factors de- 
rived from the PDR and the LVG models are very 
similar. For Henize 2-10, the filling factor is sig- 
nificantly smaller for the PDR model than for 
the LVG model. In Henize 2-10, LVG models, as 
with PDR models, do not reproduce observations 
very well due to the lack of high- J CO lines (i.e., 
12 CO(8-7) and up) which would constrain the lo- 
cation of the maximum of the CO cooling curve. 
However, compared to LVG models, PDR models 
tend to fit the series of CO lines better, particu- 
larly for NGC 253. So, we prefer the PDR values 
of the FF to those from the LVG. 

We have computed the total CO cooling from PDR 
models summing the contribution from all CO lines from 
12 CO(1-0) up to 12 CO(15-14). In the observations, we miss 
12 CO(5-4) and all CO lines from 12 CO(8-7) and up. For 
Henize 2-10, we obtained 1.5 x 10~ 8 WmT^sr^ 1 (simi- 
lar to the CO cooling from the "low Tk" LVG model 
(Tk = 50 K)). The lines which contribute the most to 
the total CO cooling are 12 CO (7-6) (17.4%) followed by 
12 CO(8-7) (15.7%) and 12 CO(6-5) (15.0%) (with the "low 
Tk" LVG model, the lines which contribute the most are 
also CO lines with J U pp er ^ 5). For NGC 253, we deduced 
a total CO cooling of 7.2 x 10 -7 Wm^sr^ 1 , (close to, 
but lower than the CO cooling from the "low Tk" LVG 
model (Tk = 70 K)). The most important line is 12 CO(6- 
5) (24.9%) followed by 12 CO(5-4) (23.2%) and 12 CO(7-6) 
(17.0%). 

For NGC 253, the cooling derived from the PDR model 
is w 40% higher than the observed CO cooling, while for 
Henize 2-10, the cooling from PDR models is a factor of 2 
larger than the measured value. This can be explained by 
the fact that, for Henize 2-10, PDR models do not fit ob- 
servations as well as for NGC 253. Actually, 12 CO(8-7) 
or 12 CO(9-8), and 13 CO(6-5) or 13 CO(7-6), would 
be very useful in order to determine the position 
of the peak of the CO cooling curve, and the opac- 
ity of the 12 CO(6-5) or 12 CO(7-6) lines. Finally, the 
strong dependence of the PDR model fit with the density 
n(H), as can be seen by comparing Fig.[3]and Fig.0 might 
also explain the difference between the two galaxies. We 
discuss this point in the next section. 

5.4. Discussion 

From the models described above, we conclude that the 
molecular gas in the starburst regions of Henize 2-10 and 
NGC 253 experience similar physical conditions: both the 
line ratios and derived properties are similar. Nevertheless, 
we can see in Fig. [S] and in Fig. El a density difference be- 
tween the two galaxies. In fact, it is difficult to fit a PDR 
model for Henize 2-10 since we do not have a CO detection 



Table 6. Parameters of the "best" PDR models. 





Henize 2-10 


NGC 253 


n(H) (cm -3 ) 


8.0±lx 10 5 


3.0±0.5x 10 b 


X (on each side) 


20000 


20000" 


v" (kms- 1 ) 


60 


187 


FF 


1.5±0.5 x 10 _a 


9.7±4 x 10-^ 



a Carral et al. (1994) found for NGC 253 a FUV flux 
of 2xl0 4 G o in units of Habing (1968) corresponding to 
1.2xl0 4 G o in units of Draine (1978); compatible with our 
value. 

b deduced from Gaussian fits to line profiles. 



Table 7. Observed and modeled line emissivity ratios for 
Henize 2-10 and NGC 253 using PDR models. 



E 


NGC 253 


HENIZE 2-10 


(ergcm~ 2 s~ 1 sr~ 1 ) 


observations 


observations 


12 CO(3-2) 
12 CO(6-5) 


0.19±0.076* 


0.22±0.083* 


12 CO(3-2) 
12 CO(7-6) 


0.27±0.11* 


0.16±0.064* 


°CO(6-5) 
12 CO(7-6) 


1.38±0.56* 


0.74±0.29* 


1Z CO(1-0) 
13 CO(1-0) 


12.73±2.5* 


16.82±6.2* 


12 CO(2-l) 
la CO(2-l) 


14.44±3.4* 


22.22±6.1* 


"CO(3-2) 
la CO(3-2) 


10.61±3.6* 


10.21±4.7* 




model 


model 


12 CO(3-2) 
12 CO(6-5) 


0.26 


0.18 


12 CO(3-2) 
12 CO(7-6) 


0.34 


0.15 


12 CO(6-5) 
12 CO(7-6) 


1.32 


0.86 


12 CO(1-0) 
13 CO(1-0) 


15.42 


15.07 


12 CO(2-l) 
la CO(2-l) 


9.72 


8.02 


12 CO(3-2) 
la CO(3-2) 


7.92 


6.27 



* : ratio derived from observations marked with asterisks in 
Table Hand Table 



for J higher than J=7-6. The reason is the large inten- 
sity of the 12 CO(7-6) line relative to 12 CO(6-5). When the 
turn over in the CO cooling curve is not well constrained, 
fitting one PDR model is not easy and the solution can be 
understood as a "lower limit" . 

Moreover lower transitions of carbon monoxide like 
12 CO(1-0) and 12 CO(2-l) were not really essential here 
since models were dedicated to the study of warm gas. So, 
the detection of high-J transitions was crucial for describ- 
ing the starburst nuclei. 

Another point is the influence of the FWHM line width. 
It is obvious for NGC 253 that the FWHM of the 12 CO(3- 
2) line (« 190fcms _1 ) is larger than that of the 12 CO(7- 
6) line (~ 95fcms _1 ), the main reason being the better 
spatial resolution of the 12 CO(7-6) data. Because NGC 
253 has a steep velocity gradient along its major axis, the 



E.Bayet et al.: The submillimeter C and CO lines in Henize 2-10 and NGC 253 



13 



convolution to lower spatial resolution, using the adopted 
source model, is expected to produce broader lines pro- 
files. We tried to take into account this phenomenon in 
our study by using LVG models with ratios in line inte- 
grated area in Kkms -1 instead of using ratios in main 
beam temperature in K. We also studied LVG models 
with a smaller Av but this implies a lower N(CO) since 
the meaningful variable in LVG models is the CO column 
density divided by the line width: N( 12 CO)/Au. So dif- 
ferent parameter pairs (N( 12 CO), Av) may reproduce the 
observations equally well. 

PDR models predict the emissivity of the CI( 3 Pi- 
3 Po) and the CI( 3 P2- 3 Pi), fine structure transi- 
tions of atomic carbon at 492 GHz and 809 GHz re- 
spectively (see Table |8"|) . For Henize 2-10, the predicted 
CI( 3 Pi- 3 Po) line emissivity is brighter than the observed 
value by a factor of 2. For NGC 253, the predicted transi- 
tions are brighter than those observed by a factor of 3 for 
CI( 3 Pi- 3 P ) and a factor of 12 for the CI( 3 P 2 - 3 Pi). We 
concluded that PDR models fitted to the CO line emis- 
sion do not reproduce the atomic carbon data very well. 
Perhaps, to better constrain models with observed atomic 
carbon transitions, we should reduce the density. The fine- 
structure lines of atomic carbon seem to share the same 
behaviour as the low-J CO transitions. This tendancy is 
found in Galactic clouds also. 

We compared observations of CO and [CI] in the center of 
the Milky Way (Fixsen et al. 1999), and in the Cloverleaf 
QSO (Barvainis et al. 1994, Tsuboi et al. 1999 and Wei/3 
et al. 2003) with our observations. These two objects 
are well known and we can compare them with 
NGC 253 and Henize 2-10, to determine the dif- 
ferences or the similarities in the physical proper- 
ties of their warm gas. We chose to plot CO lines fluxes 
(F in Wm~ 2 ) versus J U p per , instead of CO lines intensi- 
ties (I in Wm^ 2 sr^ 1 ) versus 3 U pper as in Fig. and in 
Fig- El because the Milky Way data found in Fixsen et al. 
(1999) were only available in these units. The Cloverleaf 
QSO data (in Jykms^ 1 in Barvainis et al. 1994, Tsuboi et 
al. 1999 and Wei/3 et al. 2003) have been converted to the 
same units using the Jy per K factors and the main beam 
efficiencies of telescopes (see references above). Line fluxes 
(F in Wm~ 2 ) are listed in Tableland the flux ratios in 
Table HTJ1 

CO line fluxes (F in IFm~ 2 ) for the Cloverleaf QSO, the 
nucleus of the Milky Way, NGC 253 and Henize 2-10 are 
shown in Fig.0 To compare the CO cooling curves 
more easily, the line fluxes have been scaled. The 
scaling factors are listed in the figure caption. In 
Fig. [3 PDR model predictions have been used for the 
CO lines which are lacking observations: 12 CO (5-4) and 
12 CO(8-7) for NGC 253 and for Henize 2-10 (see Table©. 
Note that the PDR model used for Henize 2-10 PDR does 
not reproduce observations as well as the PDR model used 
for NGC 253 (see Fig. Hand Fig.©. 
We compare the CO cooling curve for the center 
of the Milky Way, the Cloverleaf QSO, Henize 2- 
10 and NGC 253. The turnover positions for Henize 



Table 8. Results of PDR models for [CI] in Henize 2-10 
and NGC 253. We reported also observations from Tablc[5] 
and|3 The PDR model predictions have been multiplied 
by the line width ratio, A Av(pdr) ' an<1 by the surface 
filling factor, PDR_FF, (see Table©. 



E 


NGC 253 


HENIZE 2-10 


(ergcm~ s~ 1 sr~ ) 


observations 


observations 


CI( J Pi- J P ) 


3.5x 10- b 


3.0x 10"' 


CI( 3 P 2 - 3 Pi) 


3.2x l(T b 


not observed 




PDR model 


PDR model 


CI( J Pi- J P ) 


1.1 x 1(T 4 


5.1x 10 _v 


CI( a P 2 - a Pi) 


3.9x 1(T 4 


1.9x 10~ b 



2-10, NGC 253 and the Cloverleaf QSO are close to each 
other, near 12 CO(6-5) or above, while the turnover of the 
CO cooling curve of the nucleus of the Milky Way is near 
12 CO (4-3). We may ask if this observed difference between 
the Milky Way and the three other sources is consistent 
with a difference in linear resolution, or if it is due to dif- 
ferences in inherent physical conditions between the four 
environments. The linear resolution for NGC 253 corre- 
sponding to the adopted beam size of 21.9" is 265 pc. For 
Henize 2-10, with a distance of 6 Mpc and the same beam 
size, we obtained a linear resolution of 637 pc. For the 
Cloverleaf QSO, we computed the angular distance from 
H and q in Wei/3 et al. (2003), we obtained D A = 1056.4 
Mpc and a linear resolution of 13 kpc. For the nucleus of 
the Milky Way, the 7° COBE-FIRAS beam gives a linear 
resolution of 1035 pc in the Galactic Center (distance of 
8.5 kpc). We conclude that despite the wide range 
of linear resolutions, the turnover of the CO cool- 
ing curve is found at roughly the same J value, 
above 12 CO(6-5) for Henize 2-10, NGC 253 and the 
Cloverleaf QSO. However, at the same linear scale 
as Henize 2-10, the gas in the Milky Way nucleus 
shows less excitation. Thus, it seems likely that 
differences in the shapes of the CO cooling curves 
are due to differences in the ISM properties in the 
target galaxies. The fact that the Cloverleaf QSO, 
NGC 253 and Henize 2-10 have the same ratios is 
related to the similarity of the physical properties 
of warm gas for these three sources, which trans- 
lates into similar CO line ratios (see Table ITO|> . 

6. Conclusions: 

We observed Henize 2-10 in the rotational lines of carbon 
monoxide 12 CO(J=3-2), (J=4-3), (J=6-5) and (J=7-6) 
and NGC 253 in the 3 P2- 3 Pi fine structure transitions of 
atomic carbon [CI] at 809 GHz and in the rotational lines 
of carbon monoxide 12 CO(3-2), (J=6-5) and (J=7-6). 
We show that C cooling is less than CO cooling for 
both galaxy nuclei by a factor > 10. Among observed 
lines, those which contribute most to the CO cooling are 
12 CO(6-5) and 12 CO(7-6). 
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Table 9. Observations of the Milky Way nucleus (from Fixsen et al. 1999), NGC 253 (from this paper), Henize 2-10 
(from this paper) and the Cloverleaf QSO (from Barvainis et al. 1994, Tsuboi et al. 1999 and Wei/3 et al. 2003). 



F 

(Wm,- 2 ) 


Milky Way 


NGC 253 


HENIZE 2-10 


Cloverleaf 


CI( 3 Pi- 3 P ) 


Z T. .-. T71 

1.9x 10" iu 


4.6x 10~ ib 


„ ^ _ „ "1 y 

3.8x 10" is 


z — — tt" — m 

1.7x 10 - ^ 


CT( 3 P 2 - 3 Pi) 


— Z — T. ,, ill — 

1.9x 10" iu 


4.0x 10" lb 


not observed 


„ „ — - — <jri 

3.9x 10 Zi 


i2 CO(l-0) 


2.7x 10 _il 


2. Ox 10 -1 ' 


4.8x 10~ ly 


<1.5x 10 -24 


i2 CO(2-l) 


l.lx 10- xu 


1.9x 10 _ib 


2.5x 10~ is 




i2 CO(3-2) 


2.0x l(T iu 


4.4x 10~ ib 


9.0x 10" i8a 


4.3 x 10"^ 


i2 CO(4-3) 


3.0x 10~ iu 


1.3 x 10 _lb 


1.4x 10 _i/ 


9.11X \^'" 


i2 CO(5-4) 


2.8x 10~ lu 


2.2x 10~ ib " 


3.3x lO -17 ' 


1.3x 10~ 22 


i2 CO(6-5) 


1.9x 10~ iu 


2.2x 10~ lb 


2.9x 10 _i/ 




12 CO(7-6) 


1.7x 10 _lu 


1.7x 10~ lb 


4.0x 10~ 17 


3.6x 10~ 22 


12 CO(8-7) 


1.8x 1(T 1U 


6.7x lO" 10 " 


4.9x 10 -17 ' 





a : We chose to use 12 CO(3-2) observations from Meier et al. (2001) because they seem to agree with LVG and PDR models 
(see Fig. better than our observations. 
6 :These values are computed from PDR models. 



Table 10. Flux ratios for the center of the Milky Way, NGC 253, Henize 2-10 and the Cloverleaf QSO. 



F 

(Wm.- 2 ) 


Milky Way 


NGC 253 


HENIZE 2-10 


Cloverleaf 


CI( a P 2 - 3 Pi) 


1.00 


1.15 


not observed 


0.42 


12 CO(3-2) 
12 CO(1-0) 


7.37 


22 


18.75 


> 28.6 


°CO(3-2) 
12 CO(2-l) 


1.84 


2.32 


3.6 




12 CO(3-2) 
12 CO(4-3) 


0.67 


0.34 


0.64 


0.47 


u CO(3-2) 
12 CO(5-4) 


0.72 


0.21* 


0.27* 


0.33 


12 CO{3-2) 
12 CO(6-5) 


1.03 


0.20 


0.31 




12 CO(3-2) 
12 CO(7-6) 


1.18 


0.26 


0.23 


0.12 


°CO(3-2) 
12 CO(8-7) 


1.09 


0.66* 


0.18* 





*:These values are computed from PDR models. 



Such high-J transitions are needed to constrain the physi- 
cal conditions in starburst nuclei. We used both LVG and 
PDR models for each galaxy: the molecular gas in the 
Henize 2-10 nucleus is well described with a LVG model 
defined by AV = 60 kms" 1 , T k w 50-100 K, ^gg « 30, 
n(H 2 ) > 10 4 cm- 3 and N( 12 CO) = 3.5 ± 1 x 10 18 cm" 2 . 
For the NGC 253 nucleus, we derived the following pa- 
rameters: AV" = 190 kms" 1 , T k « 70-150 K, ^gg « 40, 
n(H 2 ) > 10 4 cm' 3 and N( 12 CO) = 1.5±0.5 x lO 1 ^ cm- 2 . 
PDR models provide equally good fits to the data. PDR 
models were used to give us more information about the 
physical parameters in these media. We succeeded in 
reproducing the observations with an accuracy of about 
20% using a PDR model defined for Henize 2-10 by 
n(H)=8.0 ± 1 x 10 5 cm' 3 and X = 20000 (we modelled 
the source as a plane-parallel slab, illuminated on both 
sides by FUV radiation). For NGC 253, we derived 
the following model: n(H)=3.0 ± 0.5 x 10 5 cm- 3 and 



X = 20000. 

Thanks to those models, we predict that, for distant 
galaxies, to obtain properties of cold and warm gas, the 
most interesting CO lines to be observed are 12 CO(1-0), 
(2-1), (3-2), (4-3), (5-4), (6-5), (7-6) and 12 CO(8-7). 
12 CO (9-8) and higher- J CO lines will be weaker and 
more difficult to detect in distant galaxies. The reasons 
are twofold. First the most intense lines are 12 CO (6-5) 
and 12 CO (7-6), second, though 12 CO (9-8) can have 
a significant contribution to the cooling, its antenna 
temperature is significantly lower than for other CO 
lines (eg 12 CO(4-3) or 12 CO(6-5)). Also, the typical T K 
and n(H2) values only allow significant excitation over 
a large fraction of the galaxy for 3 up ^10 transitions. 
For higher-J transitions the values of Tk and n(^) high 
enough to excite 3 up ^10 will be confined to a volume 
too small to produce a detectable signal. 
Data on the 12 CO (8-7) and on the 13 CO (6-5) lines would 
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Fig. 7. Flux (Wm~ 2 ) vs 
Milky Way (top left), NGC 253 (top right), Henize 2- 
10 (bottom left) and the Cloverleaf QSO (bottom right). 
PDR models have been used to obtained fluxes (plotted in 
grey) of 12 CO(5-4) and 12 CO(8-7) in Henize 2-10 and in 
NGC 253. Observations are shown in black. To facilitate 
the comparison, we applied scaling factors of 1 x 10 7 for 
the Cloverleaf QSO data, 1 x 10~ 5 for the nucleus of the 
Milky Way, and 100 for Henize 2-10. Data for NGC 253 
are shown with no scaling factor. For the Cloverleaf QSO, 
the 12 CO(1-0) line flux is an upper limit (see Tsuboi et al. 
1999). 



be extremely useful in further studies; 12 CO(8-7) will 
help to localize the maximum of the CO cooling curve, 
while 13 CO (6-5) will help constrain the warm gas column 
density. Even now, it is still difficult to detect these lines 
from the ground; with APEX receivers which may cover 
atmospheric windows up to a frequency of 1.4 THz 2 , 
these lines could become available. ALMA will give 
access to better spatial resolution, for resolving 



http:/ /www. mpifr-bonn.mpg.de/div/mm/apex. html 



individual molecular clouds : And with HIFI and 
PACS on board the Herschel satellite, lines with 
frequencies up to 5 THz will become accessible. 

We compared properties of warm gas derived from our 
models with the properties of the nucleus of the Milky 
Way and of the Cloverleaf QSO. We concluded that the 
ISM in NGC 253, Henize 2-10 and in the Cloverleaf QSO 
are similar, leading to similar CO excitation, while the 
nucleus of the Milky Way exhibits lower excitation CO 
lines. 
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